The planar pyrochlore antiferromagnet: A large- N analysis 
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We study possible quantum phases of the Heisenberg antiferromagnet on the planar pyrochlore 
lattice, also known as the checkerboard lattice or the square lattice with crossings. It is assumed 
that the exchange coupling on the square-lattice-links is not necessarily the same as those along 
the crossing links. When all the couplings are the same, this model may be regarded as a two 
dimensional analog of the pyrochlore antiferromagnet. The large- TV limit of the Sp(iV) generalized 
model is considered and the phase diagram is obtained by analyzing the fluctuation effects about 
the TV — ► oo limit. We find a topologically-ordered Z2-spin-liquid phase in a narrow region of 
the phase diagram as well as the plaquette-ordered and staggered spin-Peierls phases as possible 
quantum-disordered paramagnetic phases. 
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I. INTRODUCTION 

Recently frustrated magnets have attracted much in- 
terest of theorists and experimentalists as new kinds of 
frustrated magnetic systems have become available for 
experimental studiesi^ Armed with much progress in 
understanding the classical frustrated magnets*^ the at- 
tention is now focused on quantum counterparts. Frus- 
tration often leads to a large degeneracy of the classi- 
cal ground states and the resulting frustration-enhanced 
fluctuations can suppress classical long-range spin or- 
der, encouraging the possibility of spin-disordered ground 
states even in two or three dimensions £1 

It is not obvious, however, what kinds of quantum 
ground states may arise as a result of the fluctua- 
tions, especially for the small values of spin. Indeed 
geometrically frustrated magnets seem to allow mul- 
tifarious quantum-disordered paramagnetic phases, in- 
cluding various translational-symmetry-breaking phases 8 
and the quantum spin liquid phases with fractionalized 
excitations&2ii2iii It has been suggested that some of 
these quantum-disordered phases may also arise in un- 
derdoped cuprate superconductors, where the frustra- 
tion of the spins may occur due to the motion of doped 
holesi 2 -^ 

The possible existence of quantum-disordered phases 
on the pyrochlore lattice is a particularly interesting 
question in view of the fact that it is generally harder 
to suppress a long-range spin order in three dimen- 
sions. The progress in understanding this issue has 
been slower than the cases of two dimensional frus- 
trated magnets mainly due to the significantly bigger 
Hilbert space To circumvent this difficulty, the 
planar version of the pyrochlore lattice, also known as 
the checkerboard lattice or the square lattice with cross- 
ings, has been introducedi 17 i 18 . 19 i 20 i 21 i 22 . 23 i 24 i 25 i 26 When 
all the Heisenberg exchange couplings are the same, the 
checkerboard lattice has the same local structure of the 
pyrochlore lattice; a network of corner-sharing tetrahe- 
dra. Moreover, the size and topology of the ground state 



manifold are also identical in two cases £ Therefore one 
may expect that the studies of the checkerboard lattice 
will shed some light on the pyrochlore lattice problem. 

In this paper, we examine a Heisenberg model on 
the checkerboard lattice (illustrated in Fig^l, where the 
Heisenberg exchange couplings on the horizontal and ver- 
tical links are in general different from those along the 
diagonal links. The motivation for studying such a model 
partly comes from the fact that no symmetry of the lat- 
tice can turn the first neighbor bonds to the second neigh- 
bor bonds so that not all the bonds on its tetrahedra are 
equivalent in contrast to the pyrochlore lattice case. The 
Hamiltonian for this model can be written as 



h = Ji Sj • Sj 

(ij) 
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where Si are S =1/2 operators at the site i. Here J\ > 
and J2 > are the antiferromagnetic exchange couplings 
on the nearest-neighbor (horizontal and vertical) links 
and the second-neighbor (diagonal) links, respectively. 

We consider the generalization of the physical spin 
SU(2) = Sp(l) symmetry to Sp(N) and study the quan- 
tum ground states of the Sp(iV) Hamiltonian in the large- 
N limitiS Some of the phases obtained in the large-iV 
limit may not appear in the SU(2) model. Even these 
phases, however, may still be relevant to some physical 
systems whose microscopic Hamiltonians are "near" the 
parameter space of the original Hamiltonian. For exam- 
ple, these phases may arise when the original system is 
deformed by substitutional doping or hydrostatic pres- 
sure. 

More specifically, the first step toward the Sp(iV) gen- 
eralization starts with the bosonic representation of the 
SU(2) spin operators, S, = f&j"^^, where a, /3 =T, j 
labels two possible spin states of each boson, h ia , and 
the constraint ni, — b\ a bi a — 2S must be imposed at 
each site. Then, apart from an additive constant, the 
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FIG. 1: The checkerboard lattice. The exchange Ji acts be- 
tween sites separated by the horizontal or vertical links while 
the exchange J2 acts across the diagonal links. 



Heisenberg Hamiltonian can be written as 



H 



(2) 



where Jij = J\ on the horizontal and vertical links, and 
Jij = J2 on the diagonal links. Here e Q( a is the anti- 
symmetric tensor of SU(2). The generalization to Sp(iV) 
symmetry can be achieved by introducing N flavors of 
bosons on each site. The constraint must be modified 
to n b = b\ a b la = 2NS, where a = 1, ...,27V is a Sp(7V) 
index. For the physical case, N — 1, S takes half-integer 
values. The corresponding Sp(iV) Hamiltonian is 



H 



1 

~2N 



J2W^ a bY)(J^ 5 b i7 b jS ), (3) 



where J a ^ = J a p — —Jp a is the generalization of the e 
tensor of SU(2); it is a 2N x 2N matrix that contains N 
copies of e along its center block diagonal and vanishes 
elsewhere £ 

In the N — > 00 limit at a fixed boson density per fla- 
vor, rib/N = 2S = K, a mean field theory for S = k/2 
can be obtained and the fluctuations about the mean 
field state give rise to a gauge theory^ The mean-field 
phase diagram at N — > 00 is shown in Fig|3|as a function 
of Jij{J\ + J-i) and 1/S. At large values of 5, various 
magnetically long-range-ordered (LRO) phases appear 
and are represented by the ordering wavevector (91,92) 
measured in units of 1/ (nearest-neighbor spacing). The 
short-range-ordered (SRO) phases at small values of S 
correspond to the quantum-disordered phases with short- 
range equal-time spin correlations enhanced at the cor- 
responding wavevectors. 

One of the most notable findings in our studies is 
the existence of a quantum-spin-liquid phase with Z 2 - 
fractionalized excitations^! in a narrow region of the 
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FIG. 2: The four-site unit cell of the checkerboard lattice 
and the twelve link variables Qij. 



phase diagram (denoted by (ir,q) SRO). This region, 
Ji ^ J2 and the small values of S, has not been ex- 
plored in previous works. After taking into account fluc- 
tuation effects, the plaquette-ordered state emerges from 
the (tt, tt) SRO phase; similarly a staggered spin-Peierls 
state arises from the (tt, 0) SRO phase. The result on 
the (tt, tt) SRO phase is in accordance with the previ- 
ous finding of the plaquette-ordered phase for Ji = J2 
and small values of Si^i^^ The transition from the 
Z2 quantum-spin-liquid to the plaquette-ordered phase 
or the staggered spin-Peierls phase is described by a Z 2 
gauge theoryiSL22i22^S Our results may be also consis- 
tent with the large-.!? study of the SU(2) model with 
J 1 J 2 M> Detailed comparison with the previous works 
will be made later in this paper. 

The rest of this paper is organized as follows. In section 
II, we describe the mean-field phase diagram and various 
mean- field phases at N — > 00. In section III, the effect 
of singular fluctuations at finite TV on the paramagnetic 
phases is discussed. We explain the relation between the 
previous works and ours, and conclude in section IV. 



II. MEAN FIELD PHASE DIAGRAM 



The mean-field theory can be obtained by introduc- 
ing the directed link fields, Qij = —Qji. These fields 
are used to decouple the quartic boson interactions in S 
by a Hubbard-Stratonovitch transformation. After this 
decoupling, the effective action contains the terms 



i>3 



N\Qij\ 2 -QijJ al }b\ 



V] + c.c. 



where r is the imaginary time and the ellipses repre- 
sent standard terms which impose the canonical boson 
commutation relations and the constraint £ At the sad- 
dle point of the action, we get 



1 

N 



J U ia"jf3 
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FIG. 3: [Dashed line: first order transition. Solid line: con- 
tinuous transition.] Large-iV phase diagram of the Sp(iV) 
checkerboard lattice model as a function of J 2 /(Ji + J 2) and 
1/S. The LRO phases break spin-rotation symmetry. The 
spin order is collinear and commensurate in the (tt, tt) and 
(71", 0) LRO phases while it is helical and incommensurate in 
the (tt, q) LRO phase. The SRO phases preserve spin-rotation 
invariance. At J1/J2 = 1, we find zero expectation value of 
the diagonal bond variables at all re = 2S. The long-range 
Neel order develops above re = 2S « 0.4; below this value, 
the Neel correlations are only short-ranged. In addition, for 
large values of S, we get the Neel (n, tt) LRO for J2 < Ji and 
(tt,0) LRO for J 2 > Ji. 



The four-site unit cell of the checkerboard lattice, as de- 
picted in FigEl has twelve of these Qij fields. When S 
becomes large, the dynamics of S leads to the condensa- 
tion of the bf bosons and one obtains a nonzero value of 



(6f> 



(6) 



of link variables are nonzero and equal on the horizon- 
tal (Q) and vertical (P) links, while the values on the 
diagonal links (R) are zero. 



2. (tt, 0) and (0,7r) LRO states 

This magnetically ordered phase is characterized by 
nonzero values of (Qij) on the horizontal and diagonal 
links for the (tt, 0) LRO state or on the vertical and di- 
agonal links for the (0, tt) LRO state. There are two 
gauge-nonequivalent choices for the values of R\ , R2 , R3 
and R4. One state has R4 = R% = R3 = R4 with P = 
and Q ^ 0, and the other has R% = R$ = —R and 
R2 = R4 — R with P ^ and Q = 0. These two states 
are interchanged under a tt/2 rotation. Moreover, these 
states can be understood as diagonally coupled horizon- 
tal or vertical Heisenberg antifcrromagnctic chains. 



3. Helical (n,q) and (q,ir) LRO states 

This helically-ordered phase is characterized by 
nonzero values of (Qij) on the horizontal, vertical and di- 
agonal links. Once again, our results indicate that there 
are two gauge-nonequivalent choices for the values of R\ , 
R2, R3 and R4. One state has R\ = Ri = R3 = R4 
with Q > P and the other has R\ = R3 = —R and 
R2 = R4 = R with P > Q. These two states are inter- 
changed under a tt/2 rotation and correspond to spirals 
ordered in the horizontal or vertical directions. Notice 
that this is a long-range incommensurate spin order and 
the spin structure factor peaks at the incommensurate 
wavevector (tt, q) or (q, tt). Here q varies continuously as 
J2/J1 changes. 



Since the large- N limit of S is taken for a fixed value of 
k = rib/N = 2S, depending on the ratio J1/J2 and the 
value of k, the ground state of S at T = can either break 
the global Sp(7V) symmetry and posses magnetic LRO or 
be Sp(AT) invariant with only SRO; thus the ground state 
energy should be optimized with respect to variations in 
(Qij) and xf for different values of J2/J1 and k. We also 
find that each saddle point may be described by a purely 
real (Qij } . The resulting phase diagram is shown in Fig|3] 
We describe various magnetically ordered (xf ^ 0) and 
paramagnetic (xf = 0) phases as follows. 



A. Magnetically ordered phases 

1. Neel (tt,tt) LRO state 

This is the long-range-ordered state in which (Si) is 
collinearly polarized in opposite directions on two sublat- 
tices. In the appropriate gauge, the expectation values 



B. Paramagnetic phases 

1. Decoupled chains 

As its name suggests, this phase is a grid of decou- 
pled Heisenberg antiferromagnetic chains. In the large- N 
limit, it corresponds to the solution where the (Qij) are 
nonzero only on the diagonal links. 



2. (tt,tt) SRO state 

This state is obtained by quantum disordering the Neel 
state (it occurs at k < 0.4). The expectation values of 
Qij have the same structure as those of the Neel state. 
As described in previous works, the quantum fluctua- 
tions in this phase can be understood via a compact U(l) 
gauge theory*! This theory is always confining, thus the 
bf bosons bind to form a S = 1 quasiparticle above a spin 



4 



gap£ At finite N, gauge-theoretic consideration of singu- 
lar fluctuations lead to the plaquette order (see section 
III). 



3. (it, 0) and (0, tt) SRO state 

This is the quantum-disordered state of the (tt, 0) or 
(0, tt) LRO state and the expectation values of Qij have 
a similar structure as those of its ordered counterpart. 
Even though this state has a gap to all spin excita- 
tions, the symmetry of tt/2 rotations between the vertical 
and horizontal directions is broken. Again, the gauge- 
theoretic consideration at finite N leads to the staggered 
spin-Peierls phase (see section III). 



regime where it can be mapped onto a quantum dimer 
model, and this yields an appealing physical interpreta- 
tion of the computation. However, the mapping to the 
dimer model is not available in all cases, and here the 
instanton analysis is a bit more subtle. It turns out that 
the (tt, 0) SRO phase is one of the cases where the dimer 
model mapping is absent. Consequently, we will eschew 
the dimer interpretation here, and present details of the 
computation using the compact U(l) gauge theory. 

As has been argued elsewher o 8 ! 31 spin fluctuations in 
the paramagnetic phase of collinear antiferromagnets are 
described by the following compact U(l) lattice gauge 
theory 



z A = n 



-y^- exp JT cos ( e M^A^ iA ) 
n V □ a 



4- (ir,q) and (q,ir) SRO state 

This state is obtained by quantum disordering the 
(tt, q) or (q, tt) LRO state and the expectation values of 
Qij have a similar structure as those of its ordered coun- 
terpart. As discussed in previous works, nonzero expec- 
tation values of the diagonal link fields leads to the pres- 
ence of the "charge" -2 Higgs scalar coupled to a compact 
U(l) gauge field£i°- As a result, this phase corresponds 
to the deconfined phase of the Zi gauge theory and sup- 
ports deconfmed spinons with 5 = 1/2 above the spin 
g a p ( 27 i 28 i 29 4 30 Again, even though this state has a gap 
to all spin excitations, the symmetry of irjl rotations 
between the vertical and horizontal directions is broken. 
This phase is topologically ordered and it leads to an 
additional fourfold degeneracy of the ground state on a 
torus, 27 ! 28 ! 29 ! 30 



III. FLUCTUATION EFFECT AND 
EMERGENCE OF 
TRANSLATIONAL-SYMMETRY-BREAKING 
PHASES 

In this section, we will discuss the influence of the fluc- 
tuations about the mean-field on the spin-singlet sector 
of the (tt, tt) and (tt, 0) SRO phases. It has been known 
that the fluctuations at finite-A^ can be understood via 
the compact U(l) gauge theory^ 8 Regular perturbative 
corrections in 1/N and higher-orders do not qualitatively 
change the mean- field ground states at N — ► oo. On the 
other hand, singular effects of the hedgehog instantons 
in the gauge theory and their Berry phases lead to the 
qualitative changes, typically the emergence of various 
translational-symmetry-breaking phases^. For simplic- 
ity, we will only consider the case where 2SN is an odd 
integer. Notice that, for the physical SU(2) case with 
N = 1, it corresponds to the half-integer values of S. 

In the dual formulation, these instanton effects take the 
form of a statistical mechanical "interface" or "height" 
model. In most cases this height model has a limiting 



(7) 



Here j denotes the sites of cubic lattice in discretized 
spacetime, the subscript a the sites of the dual cubic lat- 
tice, /j,, v, A extend over x, y, r, A M is the discrete lattice 
derivative in the fj, direction, Aj^ is the compact U(l) 
gauge field, and the coupling K a can take distinct values 
of different dual lattice sites as determined by the sym- 
metry of the underlying lattice. The fixed field rjj is the 
key quantity which distinguishes different SRO phases, 
and encapsulates the staggering of the spins in the local 
collinear order. It has the values 



Vj = 



(_l)i. 



for (tt, tt) SRO 
for (tt, 0) SRO 



(8) 



Note that rjj is independent of the r co-ordinate. 

We now proceed with a duality mapping of J7J). For 
the (tt, tt) SRO phase, the procedure is essentially iden- 
tical to that in Ref. |3l], and it leads to a dual interface 
model identical to that found for the (tt, tt) SRO phase of 
the Shastry-Sutherland lattice 10 . From these results we 
obtain the plaquette-ordered state as shown in Fig. 4(a). 
Similar conclusions were also reached in Refill. 

We now apply the same procedure to the (tt, 0) SRO 
phase. The only difference here is in the value of rjj in JSJl, 
but otherwise the initial steps are the same as before 3 ^. 
First, we find an integer valued vector field, a ^ on the 
links of the dual lattice such that its curl is a vector rjj 
pointing in the r direction 



^iiv\ A„a a ^ 



(9) 



A convenient choice for a ^ is shown in Fig. [S] Then, 
we write the cosine term in Za in the Villain periodic 
Gaussian from, and perform a standard duality transfor- 
mation by the Poisson summation method. This maps 
Za to a dual interface model 



Z h = 



{h a } V a,n 



(10) 
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FIG. 4: Singular fluctuations at finite N lead to (a) the 
plaquette-ordered phase in the (n, n) SRO state and (b) the 
staggered spin-Peierls phase in the (tv, 0) SRO state. 
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FIG. 5: Specification of the non-zero values of the fixed field 
o2 M . The circles are the sites of the direct lattice, j, while 
the crosses are the sites of the dual lattice, a; the latter are 
also offset by half a lattice spacing in the direction out of the 
paper (the fi = r direction). The a^L are all zero for /j, — t,x, 
while the only non-zero values of a° ay are shown above. 



Here h a are the integer valued heights of the interface 
model. Next, we decompose into its curl and diver- 
gence free parts by writing 



(11) 



The values of the new fields ( a and y$ are shown in Fig-El 
and Fig. The values of £ will be particularly impor- 
tant to us, and these take the values (w — 1/4, (x — 
1/4, ( Y = 3/4, ( z — 3/4 on the four dual sublattices W, 
X, Y, Z shown in Fig. Inserting into i|l(J|) we 
obtain the final form of the interface model 



Z h = 




where 



H a = h a + £ a . 



(12) 



(13) 



Finally, as is usual in the mappings from interface models 
to sine-Gordon models, we can 'soften' the integer value 
constraint on the h a , and replace H a by a real valued field 
Xa which experiences a cosine potential with minima at 
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FIG. 6: The four dual sublattices upon which the height 
offsets take the values C,w = 1/4, C,x = 1/4, C,y =3/4, and 
<z = 3/4. 
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FIG. 7: Specification of the non-zero values of the fixed field 
yjfi- Only the fj, = r components are non-zero, and these are 
shown above. 



positions satisfying l|13|) . This yields finally a sine Gordon 
model for the field Xa( r ) with the action 



S x = I dr 



(ab) 



K 

I -T^Xaf ~ Va COs(27T(Xa ~ Ca)) 



(14) 



where the sum over (ab) extends over nearest neighbor 
sites, and K is the stiffness towards spatial fluctuations 
of the interface height. K T is the corresponding stiffness 
towards time-dependent fluctuations, and, for simplicity, 
its value is assumed to be independent of a. The symme- 
try of the lattice requires that the strength of the periodic 
potential take two possible values, y a — y\ or y a — y2 de- 
pending upon whether the plaquette a has diagonal Ji 
links across it or not. 

The optimal interface configurations can be deter- 
mined by the minimization of S x by a set of time- 
independent values of \w, Xx, Xy, and xz- Then, the 



6 



problem reduces to the minimization of the following en- 
ergy as a function of four real variables: 



E x = K 



(Xx - Xw) 2 + (Xw " Xy) 2 
+ (Xy - Xzf + (xz - Xx) 2 

- yi sin(27rxw0 - sin(27r\;y) 

- y 2 sin(27rxx) - sin(27rxz) 



(15) 



Our analysis here is valid only for small y\ and y 2 , so we 
will analytically determine the minima in power series in 
?/i,2- Let us define 



Xw = Xl + X2 + X3 

Xx = X1-X2+ X4 

XY = X1+X2- X3 

Xz = X1-X2- Xi- 



(16) 



At the saddle points of E x , we obtain 



\2 
X3 

Xi 



16K 2 



sin(4^xi) + 0(yl 2 ) 



^cos(2 7 rxi)+C( 2 / 1 3 )2 ) 



2K 



cos(2^xi) + 0(yl 2 ). 



(17) 



The average interface height, xi; is determined by the 
minimization of 



E x = E a + ylcos(47rxi) + £cos(87rxi) - 
where Eq is a constant independent of \l > 



(18) 
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vl) 
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3K 3 



16K 3 



(yl + yt) 
(yhl), 



(19) 



and all omitted terms are of order y\ 2 or higher. We 
now have to minimize l|18l) to determine xi- Then from 
(|17fl . we know X2,3,4) and hence the configuration of the 
interface heights. For small y\ and y 2 , A < and \A\ 3> 
\B\, therefore the saddle point solution is given by xi = 
0, 1/2. Taking X i = 0, we get X 2 = 0, X 3 = If and 
X4 = ff- Subsequently, using (JTSJ, we obtain 



Xw 



2K' 



Xx 



ny 2 
2K : 



Xy 



nyi 
2K' 



Xz 



ny 2 
2K 



(20) 



Following the previous work by Read and SachdevjSi 
the static "electric field" in the compact U(l) gauge 
theory on the links of the lattice of spins is given by 
iE a = g J2{3 £ af3A/3Xi where a,/3 = x,y. From this 
equation, we find iE x = ±gir(y 2 + y\)/(2K) depending 
on whether the electric field resides in the even (+) or 
odd (— ) row, and iE y = ±gir(y 2 — yi)/(2K) depend- 
ing on whether the electric field resides in an odd (+) 



or even (— ) column. Finally, these electric fields can be 
related to the underlying spin correlation function: this 
requires determining the analog of the couplings in (3.17) 
of Ref 32. However, the main effect can also be deduced 
here from simple symmetry considerations. The values 
of rjj in |(HJ| show that we can define an orientation to ev- 
ery horizontal link of the square lattice (from r)j = +1 to 
rjj = — 1) but not to the vertical links. Consequently, the 
uniform horizontal electric field strengthens every link 
oriented parallel to it (say) and weakens those antiparal- 
lel to. Also, there is no corresponding modulation in the 
vertical direction. Taken together, these considerations 
show that the static electric fields determined above lead 
to the staggered dimerization pattern shown in Fig. 4 (b). 



IV. CONCLUSION AND RELATION TO 
PREVIOUS WORKS 

We investigated quantum ground states of the Sp(7V) 
antifcrromagnetic Heisenberg model on the checkerboard 
lattice, where the Heisenberg exchange couplings on the 
horizontal and vertical links (Ji) are in general not the 
same as those along the diagonal links (J 2 ). The mean- 
field phase diagram of the model was obtained in the 
N — > 00 limit as a function of J 2 /(J\ + J 2 ) and l/S. 
We found various magnetically-ordered phases at large 
values of S, and the corresponding quantum-disordered 
paramagnetic phases at small values of S. The effect of 
singular fluctuations about the mean-field states was ex- 
amined in the quantum-disordered paramagnetic phases 
and it leads to the emergence of various translational- 
symmetry breaking phases. 

At small values of S, when J 2 is small or compa- 
rable to Ji, the ground state is the plaquette-ordered 
phase (shown in Fig. 4 (a)) that is obtained by quantum- 
disordering the Neel ordered state; this is as predicted in 
a related study on a lattice with a similar symmetryifi, 
and agrees with the studies of the J\ = J 2 case by 
the exact diagonalization and large- N studies ^i^S^ As 
J 2 is increased, a topologically ordered Z2-spin-liquid 
phased appears in a narrow region of the phase dia- 
gram. This phase is the quantum-disordered phase of the 
(71-, q) LRO phase. If J 2 is further increased, the ground 
state becomes the staggered spin-Peierls phase (shown in 
Fig. 4 (b)) that is obtained by quantum-disordering the 
(ir, 0) LRO phase. The transition from the plaquette- 
ordered (or the staggered spin-Peierls) phase to the Z 2 - 
spin-liquid phase should be described by a Z 2 gauge 
theoryi 8 i 10 i 19 i 27 i 28 i 29 i 30 On the other hand, in the J 2 > Ji 
limit, the mean-field ground state is the decoupled-chains 
SRO phase. At finite N, a weak coupling between the de- 
coupled chains will be generated and this may lead to the 
"sliding-Luttinger-liquid" phase previously discovered in 
the study of the J\ ^ J 2 case with S = 1/2; Our large- iV 
study cannot capture this physics. 

When S is very large, there exist two magnetically or- 
dered phases; (tt, it) LRO for J 2 < J\ and (71", 0) LRO for 
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J2 > Ji . This result is consistent with the recent large-5 
study of the Ji ^ Ji case^ At intermediate values of 
S, however, there exists an additional magnetically or- 
dered phase, the (tt, q) LRO phase, between the (tt, tt) 
and (tt, 0) LRO phases. This incommensurate ordered 
phase is the parent state of the iv2-spin-liquid phase. 

Finally, it is intriguing to see the emergence of the spin- 
liquid-phase in the planar pyrochlore (albeit in a narrow 
region of the phase diagram) in view of the fact that the 
existence of a three dimensional spin-liquid-phase was 



recently proposed in the easy-axis limit of the antifer- 
romagnetic Heisenberg model on the three dimensional 
pyrochlore latticed 
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